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EXTENSION OF LEADING -EDGE -SUCTION ANALOGY TO WINGS 
WITH SEPARATED FLOW AROUND THE SIDE EDGES 
AT SUBSONIC SPEEDS 

By John E, Lamar 
Langley Research Center 

SUMMARY 

This paper presents a method of determining the lift, drag, and pitching moment for 
wings which have reattached flow aroimd the leading and side edges. The method is an 
extension of the leading-edge -suction analogy of Polhamus applied to the side edges. The 
value of the term associated with the lift contribution of the reattached flow from the side 
edges Ky^gg has been found to exceed that of the leading edge for rectangular wings of 
aspect ratio less than 2. Limiting values of have been determined for rectangular 

wings. Comparisons of the results of this method with experiment indicate reasonably 
good correlation of the lift, drag, and pitching moment for a wide planform range. The 
agreement of the method with experiment was as good as, or better than, that obtained by 
other methods. 


INTRODUCTION 

Many current and proposed aircraft and missiles designed for high-speed flight 
employ highly sweptback and tapered low -aspect -ratio wings with sharp or thin edges. 
These planforms exhibit flow separation along the leading and side edges followed by 
subsequent reattachment downstream or inboard, respectively, over a large angle -of- 
attack and Mach number range. However, the effect of this separated-flow (commonly 
termed vortex flow) phenomenon is more important at subsonic speeds because of its 
larger contribution to the total aerodynamic characteristics. Polhamus in references 1, 
2, and 3 has shown that, for a variety of pointed-tips planforms and Mach numbers, the 
contribution of the leading -edge vortex to the lift and drag can be accounted for by what 
is termed the leading-edge -suction analogy. The application of this analogy requires 
only that the attached-flow leading-edge suction, available from inviscid theory, be known 
accurately. 

It is evident that the separated flow around the tips or side edges of swept and 
unswept lifting surfaces has essentially the same behavior as that around the highly 
swept leading edge. Consequently, it should be possible to predict the effect of side -edge 



vortices on the aerodynamic characteristics of low -aspect -ratio planforms by an analy- 
sis similar to those used to treat the leading edge. The purpose of the present paper 
will be to detail how one such extension to the leading -edge -suction analogy can be 
effected. The one new element required in this extension is the attached-flow side -edge 
suction force and its derivation is presented. 

Total -force and total -moment predictions, including both the leading- and side -edge 
suction-force contributions, have been made for a variety of planform shapes and are com- 
pared with previously published data and new subsonic wind-tunnel data. Also, a discus- 
sion of some other methods which attempt such predictions (refs. 4 to 11) will be given 
and comparisons with the present technique are made where possible. 

SYMBOLS 


Values are given in both SI Units and U.S. Customary Units. The measurements 
and calculations were made in U.S. Customary Units. 


AR 


B](Xi) 


b 

Cd 


AC 


aspect ratio 

coefficient of at xj of the spanwise curve fit of equation (12) 


span, cm (in.) 


drag coefficient. 


Drag 

•looSref 


Cp) 0 experimental value of drag coefficient at Cj^ = 0 


D 


_ 9Cj^ 
■^a da 


drag coefficient due only to lift, Cj) - Cp^p 
Lift 


lift coefficient. 


QooSref 


'm 


pitching -moment coefficient about the reference point, unless otherwise 

stated 11 is located at Pi tching moment 

4 qooSrefCref 




normal-force coefficient, l^ X gial force 

QooSref 
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c 


streamwise chord, cm (in.) 


c(ri) streamwise half chord at r], cm (in.) 


d(??) x-location of local midchord with respect to half root chord, cm (in.) 

F tip suction force from one side edge, N (Ibf) 

G(x) defined by equation (2), m^/^^ec (ft^^^/sec) 

j index ranging from 1 to p 


Kp = 


Q(^N,p) 

8(sin a cos a) 


^v,le " 


2 Leadir^-edge suction force from one side 

goo^ref 

9 sin^ a 


K 


■v,se 


2 Tip suction force from one side edge 

qopSref 

9 sin2 0! 


Kv,tot ^^v,le ^v,se 
M Mach number of free stream 

N maximum number of chord loadings in modified Multhopp solution 


n 


Qoo 


index ranging from 0 to N - 1 

coefficient of chordal loading function in modified Multhopp solution, m 


q^o free -stream d5mamic pressure, N/m^ (Ibf/ft^y 

R Reynolds number 

Sref reference area, m^ (ft^) 



s(x) distributed edge suction force (see eq. (1)), N/m (Ibf/ft) 

U free -stream velocity, m/sec (ft/sec) 

u induced velocity in the X-direction at a point (x,y), m/sec (ft/sec) 

V induced velocity in the Y-direction of a point (x,y), m/sec (ft/sec) 

x,y distances from a coordinate origin located at the leading -edge apex; 

X positive downstream and y positive toward right wing tip 

Ax distance along tip chord, cm (in.) 

a angle of attack, degrees 

13 = /l - m2 


y 

6 

V 

9 

0Xi 

A 


accumulated circulation at a point (x,y) (see eq. (5)), mV sec ^ftVsecj 
distributed bound vorticity at a point (x,y) (see eq. (6)), m/sec (ft/sec) 
distributed trailing vorticity at a point (x,y) (see eq. (4)), m/sec (ft/sec) 
nondimensional spanwise variable, 2y/b 
chordwise angular variable (see eq. (9)), degrees 
9 value which yields xj in equation (9) 
leading-edge sweep angle, positive for sweepback, degrees 
taper ratio, ct/Cr 


P 

a 


fraction of local chord 
density, kg/m^ (slugs /ft^j 

trailing -edge sweep angle, positive for sweepback, degrees 
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Subscripts; 


c centroid 

i particular item of location 

le leading edge 

n notch 

p potential or attached flow 

r root 

ref reference 

se side edge 

t tip 

te trailing edge 

tot total 

vie vortex effect at the leading edge 

vse vortex effect at the side edge 

THEORETICAL DEVELOPMENT 

The concept embodied in the leading-edge -suction analogy of Polhamus is developed 
in reference 1. However, to aid in illustrating the application to other edges, the pri- 
mary ideas are briefly reviewed. 

Wings which have attached flows develop suction forces along their leading edges 
if the stagnation surface does not lie along that edge. This suction force can be envi- 
sioned as arising by either of two processes: (1) the pressure near the leading edge 
acting over the edge thickness or (2) the product of the square of the induced tangential 
velocity and the distance to the edge. For a wing of infinitesimal thickness the induced 
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tangential velocity approaches an infinite value of u(x,y) as shown in sketch (a); 
however, its product (described above) is still finite. 



SIDE- EDGE 
SUCTION FORCE 


SIDE- EDGE 
SUCTION FORCE 

Sketch (a) 


If the flow separates from the wing in going around the leading edge due to its 
sharpness or thinness, or due to a combination of thickness and angle of attack, the suc- 
tion force in the chord plane is lost. However, if this separated flow forms into a shed 
vortex which causes the flow to reattach to the leeward surface of the wing, then the 
energy redistributes on the upper surface near the leading edge and consequently the 
force acts in the normal-force direction. By making the edge sharp this additional nor- 
mal force can be generated at almost all angles of attack. 

According to the analogy, the reattached line or details of the pressure field need 
not be known in advance in order to determine the reattached-flow force. However, if 
pitching -moment estimates are needed the distribution of the reattached force must be 
known. The centroid of the leading -edge suction has been used as the longitudinal loca- 
tion of this force. This assumption does not have provision for angle -of -attack effects 
on the location of the reattachment line or vortex core, hence the core is assumed 
to remain stationary near the wing leading edge. 

From the above outlined ideas it can be seen that the conditions which lead to this 
additional normal force would not necessarily be limited to wings with separated flows 
around the leading edge but could be applied to any similar situation where, in potential 
flow, an edge suction force would be produced. Sketch (a) also shows, for example, that 
along the side edge of a finite streamwise tip chord large values of v(x,y) are produced 
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due to flow around the side edges. These in combination with the infinitesimal thickness 
that lead to them would produce a finite suction force in the Y-direction (side -edge suc- 
tion force, plus or minus depending on the edge). Hence, all that is required to employ 
a generalized suction analogy to the side edges is to determine the attached-flow side- 
edge suction force. 

A mathematical procedure for computing this side force is developed in the follow- 
ing steps. Figure 1 illustrates graphically selected steps for a t5q3ical wing at a partic- 
ular chordwise and spanwise location. 

(1) The suction distribution along an edge per unit length is obtained from refer- 
ence 12 to be of the form 

s(x) = P7 tG(x)2 ( 1 ) 

where for side -edge or tip suction the term G(x) is interpreted as 

G(x) = \/| lim Vl - V v(x,y) (2) 

V l /j— 1 

where v(x,y) is the perturbation velocity in the Y-direction. 

(2) The velocity v(x,y) is related to the trailing vorticity by 

v(x,y) = I 6(x,y) (3) 

(3) The trailing vorticity of a particular x -position x^ is determined by 

(4) 

^ ^ ' 8y b ar? 


where 


pXi pXi 

r(xi,y) = J y(x,y)dx = 2 j u(x,y)dx 
^le ^le 


(5) 


See figure 1. 

(4) The subsonic bound -vorticity distribution is represented herein by 

N-1 

_ U 1 2 Y cos ng + cos(n + 1)0 
2 clri) ir L sin 9 


n=0 


( 6 ) 


which is employed in a new version of the modified Multhopp lifting surface solution of 
reference 13. 
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(5) Upon obtaining solutions for the 


equation (5) becomes 




qn(^) 

qoo 


terms from the method of reference 13, 


r^i N-l 


U TTc{rj) 


xie n=0 


y cos ng + cos(n + 1)0 qn(^) ^ 
Z; sin 6 q^ 


(7) 


or 


£I^kZl= 1 

U 7T 



^ [cos n0 + cos(n + 1)^ (jg 


( 8 ) 


where 


and 


Cy, 

X = - c{rj) cos 9 + d(? 7 ) 


dx = c(/j) sin 9 d0 

Upon integrating, the result of equation (8) is 


(9) 


( 10 ) 


_ 1> 

) y qn(^) 

sin n0 

U 7T^ 

\A qoo 

n 


n + 1 


^Xi qo(^) 

0 


qOWrn • ^ 

+ » + sin 0 

qoo 



( 11 ) 


(6) Knowledge of the r(xj,y) distribution at discrete points is not sufficient for the 
present analysis since it is the variation which is needed. This requires a 

^ Q iv) 

knowledge of the continuous variation of the — terms. In the Multhopp solution the 

Qoo 

spanwise distributions of these terms are assumed to be composed of a sine series, 

r(x* 77 ) 

which is expressible as an even power series in 77 . Hence, the ^ values com- 
puted by equation (11) will be curve fitted at each xj location with 

Bl(xi) + B2(xi)r;2 + . . . + Bp(xi)?72(p-1) (12) 




where the Bife') values are determined in the fitting process. Four terms in the series 

' rfx’ 77 ) 

were determined to represent adequately the ^ distribution in a least -squared 

sense. 
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( 7 ) Hence, the differentiation of equation ( 12 ) with respect to 77 leads to 


= I u|^/rr;?|2B2(xi)77 + 4B3(xi)?73 + 6B4(xi)77' 


\/l - 772 


Bi(xi) + B2(xi)772 + B3(xi]7?^ + B4(xi)7? 




( 13 ) 


As 77—1 the first group of terms — 0 and the second group — -00, Therefore, it is 
only these last terms in combination with their multiplier which can contribute to the 
suction. Setting 77 = 1 in all the terms except those multiplied by ^ results in 


llm = llm I 


7-1 77-1 I 2V1 - 


Bl(xi) + B2(xi)T72 + B3 (xj) 774 + B4(xi)776|> (14) 


(8) Substituting equation ( 14 ) into equations ( 4 ) and ( 2 ) results in 


= -V2 2 l(^i) + ®2 W + B3(xi) + B4(xi) 


= ..IhRUrr, "J? 

7 ]‘ 

Squaring and simplifying leads to 

t2 


G(xi)2 = mi 
16 




( 9 ) Combinii^ equations ( 1 ) and ( 16 ) yields 

I ~l2 


■N _ pU^bTT 


5 (xi) = 


16 


j=l 


or 


( 15 ) 


( 16 ) 


( 17 ) 


_ qoo^ 


i(Xi) = 


8 


1 

j=l 


( 18 ) 
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(10) The magnitude of the tip suction force from one edge can be found by 


F = r*®’* 

Jx s(xi) dx 

^le,t ^ ^ 


(19) 


or, in coefficient form 


rx, 


bTT 


^loo^ref ^^ref 




te,t 


Te,t 


1 

i=i 


dx 


(20) 


(11) Inasmuch as s(xj] depends on [y(x,y)] and since v(x,y) is linearly 
dependent on the values of which are dependent on sin a or a (in radians) for 


Q.OO 


small angles of attack, then s(x^) depends on sin^ a. Hence, the partial derivative 


with respect to sin^ a can be taken of — E — to determine a Kv-t3T>e term analogous 

qooSref 

to that obtained for the leading-edge -suction force in reference 1. Therefore the equa- 
tion defining Ky se written 


K 


2F 


v,se 


ISco^ref] 


d sin^ 01 


(21) 


With the determination of the attached -flow side force and by use of the generalized 
suction analogy, the magnitude of the additional normal force associated with the sepa- 
rated flow around the side edge can be found and its contribution to the static longitudinal 
aerodynamic characteristics can be taken into account. For example: 


Cl = Kp sin a cos^ of + (Ky^^g + ^v,se)®^*^^ cos a 


( 22 ) 


or 


Cl = Kp sin a cos^ a + Ky tot ck cos a (23) 

ACp) = Cl tan o; = Kp sin2 a cos a + Ky^tot ^ (24) 

and 

Cm = Kp sin o? cos a ^ + Ky^ig sin2 a + Ky^gg sin2 q- ^ (25) 
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where the particular x-terms equal 
Xref ■ ^c,i 

Appropriate simplifications can be made if only the leading or side edges are sharp. 

The procedure is computerized and is available from COSMIC as NASA Langley 
computer program AOS 13. 


BEHAVIOR OF Ky^se 

The behavior of Ky^se studied by examining the important parameters 

which affect its variation. They must include planform variables along with Mach number 
since both affect the span loading and, consequently, Ky gg. The planform variables 
associated with unswept or swept wings having streamwise side edges will be considered. 

Unswept Wings 

Straight trailing edge. - The geometric parameter defining rectangular wings is, of 
course, the aspect ratio. Values of K were computed based on the described pro- 
cedure and the results are presented in figure 2 as a function of this parameter. The fig- 
ure shows a rapid decrease of (solid curve) with increasing /3AI to the extent 

that for jSAl = 2 the value of g^ was only approximately one-half of the value at 
/3AI = 0. The procedure by which the limiting values of K are determined is pre- 

V j o t? 

sented next. 

The K , limit as Al — « is easy to establish because with the reference area 
v,se 

also — and the side edge forces remaining finite, K — 0. However, in order to 

V ^ O C 

establish the g^ limit as Al — 0 it is necessary to refer to reference 14 where 
certain stipulations are given for wings of very low aspect ratio. They are: (1) no chord- 
wise loading past the point of maximum span and (2) the span loadings are elliptical. For 
a wing with an unswept leading edge and streamwise tips the maximum span occurs at the 
leading edge which according to the first stipulation would require the chord loading to be 
impulsive. With the form of the loadings now specified a procedure similar to that already 
detailed can be used to determine this K value. By employing the Ct expression 

V j o -l-J 

'TT 1 

consistent with the prescribed loadings, , then se ~ ^ 

Sion of Helmbold (ref. 15) which is 


^It should be noted that for Al = 0 this is the same value as for the delta wing 
Ky le given in reference 1. That they are the same is not surprising since the delta 
wing has its sweep angle -90° as Al - 0. Hence, the leading edges are becoming like 
side edges and with Ky jg - 0 for unswept wings, their totals are identical. 
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■ 


2n/Ra 


i 


+ A + 2 


was used instead, then the values of Ky^se weighted by an expression which is 
valid over a wider aspect-ratio range. The results of these computations peld the 
long dashed curve in figure 2 which in comparison with the solid curve shows mostly the 
effect of going from an impulsive -tj^e chord loading to a finite -aspect-ratio-wing chord 
loading. 

To demonstrate the reduction in Ky^se which occurs in going from an impulsive 
chord loading to one that is either of the cot ^ type or constant (see fig. 3), while keep- 
ing the span loading elliptical and the aspect ratio near zero, computations were made 
using the same procedure as before and resulted in 

^v,se “ 1.93 

for the cot ^ loading, and 
^v,se “ 1*33 

for the constant chord loading. 

Figure 4 shows some of the results of a numerical experiment designed to deter- 
mine the aspect -ratio limit of the present method as iR — 0. From this experiment 
there are several important points to be made in this regard which follow: 

(1) The three-dimensional span loadings of wings in this aspect-ratio range are 
found numerically to be nearly elliptical, hence the effects seen here are primarily 
associated with the change in chord loadings. 

(2) The validity of the tip-suction distributions associated with the iR = 0,05 and 
iR = 0.0001 wings is suspect since these wings are outside the aspect -ratio range where 
the attached-flow method employed would yield reliable results. The validity is ques- 
tioned because the near - and far -field drag results did not agree well indicating an 
inability of the method to get good chord loadings at these jR values, 

(3) The iR = 0,1, 0,2, 0,3, and 0.5 wings had acceptable agreement between the 
two induced drag results, hence the chord -loading solutions can be assumed to be theo- 
retically correct. Because of the above points the present theoretical prediction of the 

K shown in figure 2 has been adiusted between /3iR = 0.10 and 0 to fair into u 
V,SG 

at ^jR = 0 by extrapolating the theoretical curve to the /3R = 0 value from slender - 
wing theory. 

Figure 5 presents the variation of j3K^, /3K^ and with /3R. The 

/3Kp curve is well known ^as and will not be discussed. It is interesting to 

note that at )3R = 2 the value of /3K , is about equal to the computed value of K 

v,ie ^ v,se 
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in figure 2, and with increasing l3^, the value of /BKy^le becomes larger than 
The total of /SKy^le gg stays within a ±10-percent band about tt for the 

l3M range considered; consequently, to a good approximation, g© " le 

this ^JR range. 

Notched trailing edge . - The discussion so far about the general behavior of )3Kv ge 
with /3^ and its limits for rectangular wings is applicable to wings with notched trailing 
edges. Some variations occur in between the limits reflecting the changes in side- and 
trailing-edge geometry. A reason for considering wings of this t 5 q)e was to determine if 
more side -edge vortex lift could be obtained for a given span and area. 


The major difference in the computation of the tip suction for these wings, as com- 
pared with those having straight trailing edges (although different, it poses no problem), 
is that these wings can accumulate all of the chord load at the inboard span stations (due 
to shortened chord) ahead of the tip trailing edge. This makes r(x,7j) constant for all 
chord stations aft of the notch in the spanwise curve fitting of r(x,7j) - a part of the tip- 
suction computational process - and is illustrated in figure 6. ^Because the potential- 
flow solution uses only a continuous chord description, the inner streamwise edge of the 
notch is not treated streamwise but at some sweep angle. This causes, in figure 6, 
the *^ot to be totally unchanged aft and inboard of the notch.^ 

One aspect of the notched -wing edge suction that tends to be confusing is the lack 
of edge suction along the inner edge. To understand this lack of suction it must be 
remembered that for a suction force to be present along an infinitesimally thin edge the 


3r(x.,r]) 

tangential flow at the edge must approach an infinite value. This means that 

dr] 


/ , 9r(x.,77) 

or V fXp77 j — ±°o along the inner edge, but figure 6 clearly shows that ^ only 

— -°° as T] — 1. Hence, no inner edge suction is produced for the notched ^ng 

examined. 


In order to show an example of the changes which can occur in Ky gg and Ky tot 
as a result of forming a notched wing (see sketch (b)) from a rectangular one by extending 
the side edge 10 percent and removing area at the trailing edge to keep the same aspect 
ratio the following table was prepared for M = 0: 


Planform: Al = 1 

^v,se 

^v,tot 

Rectangular 

2.1255 

2.9071 

Notched 

2.3863 

3.1650 

Percent change 

12.3 

8.9 
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Sketch (b).- Rectangular and 
notched wings. ^ = 1. 


For this aspect ratio Kp and Ky jg slightly in going from a rectangular to a 

notched wing, however, for a > 6° the larger Ky^tot the notched wing is sufficient 
to cause its Cl value to exceed that of the rectangular wing and thereafter remain 
larger. 


Swept Wings^ 

The Ky se variations for wings with leading -edge sweepback, stream wise tip 
chords, and unswept or swept trailing edges called cropped deltas and either cropped 
diamonds or arrows, respectively, are considered. 

Cropped deltas . - Cropped-delta wings are defined in terms of aspect ratio, leading- 
edge sweep, and taper ratio. The effect of varying these geometric parameters and 
Mach number on Ky^se Ky^le is given in figure 7 and their sum at M = 0 

appears in figure 8. From figure 7 the following effects were noted for three separate 
sets of independent variables: 


Set 

X 

A 

JR 

M 

^,se 

^v,le 

1 

Constant 

Increasing 

Decreasing 

Constant 

Increasing 

Decreasing 

2 

Increasing 

Constant 

Decreasing 

Constant 

Increasing 

Decreasing 

3 

Constant 

Constant 

Constant 

Increasing 




Although there is a wide variation of Ky^ge and Ky^ig the geometric parameters, 

shown in figure 7, their sums fall within a 10 -percent band around n (see fig. 8). 


^Reference 16 found two separate vortex systems for a rectangular wing, but for 
swept cropped wings the tip-vortex system merges with that from the leading edge to 
form only one system at a low angle of attack. 
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It is interesting to note that the number n keeps appearing in suction analogy 
considerations - as a limiting value at ill = 0 (Ky gg for unswept -leading -edge wings 
and Kyje delta wings), as an approximate total of Ky^se and Ky^le the 
unswept -leading -edge wings and cropped deltas over a range of aspect ratio, or as half 
the two-dimensional value of Ky^le. 

Cropped diamonds and arrows . - The cropped-diamond or cropped-arrow wings are 
defined by the same geometric parameters as the cropped deltas with the addition of 
trailing-edge sweep. In order to isolate the trailing -edge -sweep effect on Ky^se the 
following study was performed. The wing leading -edge sweep and tip-chord length were 
fixed for four different combinations of these parameters and the trailing-edge sweep 
angle was generally varied from 60° to -60°. The actual tip chord was held con- 
stant rather than \ since Ky^se has been shown in the section on unswept leading 
edges and notched trailing edges to be almost directly related to changes in tip-chord 
length. Furthermore, with varying, the taper ratio defined in terms of Cy would 
be varying for a fixed c^-. Hence, the taper ratio is only used at = 0° to fix the c^. 
The total of Ky gg and Ky j^g is presented in figure 9 rather than the separate coeffi- 
cients since, along any one of the three curves with a nonzero tip chord, both Ky ge and 
Ky,ie have the same general behavior with aspect ratio. This is because the reference 
area used to determine the coefficients is the true area which varies widely over the 
-range. In fact with the leading -edge sweep and tip chord fixed as — -90° the ref- 
erence area — « which causes Ky^tot ^ iR — 0. Thus, any finite tip or leading - 
edge-suction force will yield a zero coefficient. 

The variations of Ky tot with three separate sets of geometric parameters are 
given in the following table: 


Set 

m 

A 

°t 


liM 

n 

Increasing 

Constant 

Constant 

Increasing 

Increasing 

2 

Constant 

Constant 

Increasing 

Increasing 

Increasing 

3 

Constant 

Increasing 

Increasing^ 

Increasing 

Increasing 


^Although Ct is increasing, X is constant for unswept trailing 


edge. 

As an exercise Ky tot values were recomputed for wings with A = 63° and 
^fi=0 ~ over the -range using as a reference area that of the = 0° wing. The 
resultant variation of Ky^tot was less than 10 percent. 
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MODELS AND TEST CONDITIONS 


In order to provide needed data for comparison with the present method, an exper- 
imental investigation was conducted. A total of 10 thin flat wings were tested, each 
having stream wise tips and symmetrical beveling along the leading and side edges. Six 
of the wings had unswept leading edges along with symmetrically beveled trailing edges. 
Geometrical descriptions of the models are given in table I. Figures 10 and 11 are 
photographs of the models and typical model installations are shown in figures 12 
and 13. 

The models were tested in the Langley high-speed 7- by 10-foot tunnel at M » 0.20, 
^ 2394 N/m^ (50 psf), and R ~ 4.265 x 10®/m (l.30 x loVft). Number 80 transition 
grit at a density of ~1 particle/2 mm was added to both sides of each model approximately 
2.54 cm (1 in.) behind the wing leading edges and extending to the tips. A bolt-on bal- 
ance housing was mounted on one side of the model with a dummy housing mounted sym- 
metrically on the other side. (See figs. 12(b), 13(a), and 13(b).) The drag data have 
been corrected to a condition of free -stream static pressure acting on the base of the 
balance and dummy housings and in the chamber. No blockage or wall -interference 
corrections were made due to the small models employed at a low test Mach number. 

EFFECT OF EDGE SHAPING ON RECTANGULAR-WING RESULTS 

The results predicted by the generalized suction analogy (present method) will be 
compared with experimental data for the various models of the investigation. However, 
before these comparisons are presented the results of a limited study concerned with 
the agreement which can be obtained with wings having different edge conditions will be 
discussed. Also, the effect of separation occurring primarily at the side edges and not 
the leading edges will also be discussed. 

Figure 14 shows the effect of edge shaping on the variation of Cj^ with a for an 
Al = 2 rectangular wing. There is not much effect noted (probably due to the low test 
Reynolds numbers) except for the unsymmetrical beveled edges which exhibit a camber 
contribution to Cj^. (Some of the data were converted to Cj^ from Cl values at 
zero leading-edge thrust.) The best agreement between the present method and previous 
experimental data is achieved with the thin rounded-plate model until the leading-edge 
flow no longer reattaches to the wing. Reference 16 considers this to occur at a > 15°. 
The symmetrically beveled wing achieved the highest test Cj>^ values and was well 
predicted even at of > 15^. 

Reference 17 shows that for an Al = 0.2 rectangular wing the higher values of 
Cl are obtained for a < 37° with a thin flat plate. Hence the best agreement with the 
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present method is expected when applied to thin flat wings with symmetrical beveling 
along the leading and side edges, (It should be noted that the models constructed for the 
purpose of providing data for comparison with the present method had these features.) 

Wings with well-rounded leading edges but having separation vortices along the 
side edges (due either to symmetrical beveling or attached-flow breakdown due to 
exceeding a critical Reynolds number) provide interesting data for comparison with the 
present method. This is because all nonlinear effects, which tend to increase the curve 
of Cl plotted against a at a rate faster than that for the attached flow, must be 
attributed to the side -edge suction for angles of attack below leading -edge separation. 
Such data was found in reference 6 for rectangular wings of various aspect ratio and 
compared in figure 15 with the present method. In general, the agreement is good until 
the leading -edge flow separates and consequently affects the side -edge vortex effects. 

Not only does edge shaping affect and Cl but the drag coefficient due only 

to lift ACj) as well. Figure 16 shows the variation of ACp with aspect ratio for 
different edge conditions at four different values of Cl- As expected, for > 1.0 all 
the curves for and attached flow are slightly lower than those for wings with 

round leading edges and sharp side edges. They also show that over the entire ifl and 
Cl ranges values of ACj) due to the round-leading -edge, sharp -side -edge configura- 
tions, which retain their leading -edge suction, are lower than those for the configurations 
with sharp leading and side edges. In general, at any Cl for < 1 the value of 
ACj) for the round -leading -edge, sharp-side -edge configuration is lower than either 
attached -flow drag value due to the nonlinear contribution of Cl from the tip. Fur- 
thermore, for Cl = 0.05, 0.20, and 0.35 there are aspect ratios for which the values of 
ACj) for the configurations with sharp leading and side edges are lower than those for 
attached flow, reflecting again the increasing contribution of the side edges as Jfi 
decreases. 

COMPARISONS WITH EXPERIMENTS 
Unswept Wings 

Straight trailing edge .- Figure 17 presents some static longitudinal aerodynamic 
characteristics plotted against Cl for five rectangular wings ranging in aspect ratio 
from 0,2 to 3.0 at a Mach number of 0.2. 

In figures 17 to 21 the curves of Cq for zero percent leading -edge suction are 
computed by C^ = C^ q ® tot 100 -percent-suction 

curves are determined by Cjj q + For these five models some observations 

can be made. They are: 
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(1) The predicted Cl variation with a is in reasonable agreement with the 
experimental data. If the tip -suction contribution to Cl had not been included in the 
theory, the agreement would have been poor except at /R = 3. 

(2) The Cj) data are well predicted by the zero percent leading -edge suction and 
show the reduction m Cpj below that for attached flow at iR = 0.2, 0.3, and 0.4, which 
is consistent with the results of figure 16. 

(3) The Cm trends are reasonably well predicted for the R = 0.2, 0.3, 0.4, and 
1.0 wings although the predicted values for all except R = 1.0 show more nose -down 
moment than measured. However, it can also be seen that if it were not for the C^^ 
contribution from the side edge, the predicted moment would go the wrong way. 

(4) At R = 3.0 both the Cl variation with or and Cm data show that for 
Cl > 0.4 the leading -edge vortex lift is being lost and also, of course, its contribution 
to Cm* This occurs, as for the R = 2 wing discussed earlier, because the leading- 
edge flow no longer reattaches to the wing for a ~ 8°. With the loss in Cj^ the 
Cm data now show primarily the effect of the side-edge separated flow. 

Notched trailing edge . - The data (again static longitudinal aerodynamic character- 
istics plotted against Cl) for th® notched rectangular wing presented in figure 18 are 
predicted well by the present method and illustrate again the significant contribution that 
tip suction makes for wings in this aspect-ratio range. In comparing the data for the 
R = 1 rectangular wing with these data it can be seen that (as discussed in the section 
on the behavior of Ky gg) the values of Cl are generally higher for the notched wing 
although the lift distributions produce approximately the same Cm variation. The pre- 
dicted Cm for the notched wing is more negative than that for the rectangular wing 
(see fig. 17(d)) and only slightly overpredicts the experimental nose -down moment. 

Swept Cropped Wings 

Taper -ratio variation (deltas) . - Figure 19 presents static longitudinal aerodynamic 
characteristics as a function of Cl for a series of cropped delta wings with variable 
\ at M = 0.6 and A = 63°. The experimental data are for models which did not have 
sharp edges (ref. 18). However, the comparison shows generally good agreement 
between the data and the present method. The airfoil employed was the NACA 63A002 
which at low Reynolds number develops separated flow due to the small leading -edge 
radius. At the higher angles of attack vortex breakdown must be occurring ahead of the 
trailing edge to account for the decrease in Cl and increase in Cp. The curves of 
the various contributions of Cl plotted against a show that with increasing X the 
Cl vse contribution increases, as it should. The underprediction of Cl at the higher 
angles of attack and taper ratios, before the occurrence of vortex breakdown, can be 
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attributed to the actual shed-vortex system which produces additional induced effects 
above that predicted for an assumed core lying in the vicinity of the leading and side 
edges. An underprediction of Cl was also found for sharp-edged diamond wings in 
reference 19. 

For delta wings of Al s 2, the centroids of leading -edge suction and attached- 
flow lift are approximately the same (ref. 20); hence using the attached-flow lift centroid 
to concentrate all of the lift for the moment prediction was an acceptable procedure. 

This is no longer true if the delta-type wing has a cropped tip for then each individual 
loading centroid must be used. 

Trailing -edge sweep variation (diamond and arrow) .- Figure 20 presents some 
static longitudinal aerodynamic characteristics plotted against Cl for a family of 
cropped wings all having A = 63° and the same tip -chord length and tested at M = 0.2. 
The models differ in their f2-values which vary from -40° to 40°. For = 0° and -40° 
Cl is under predicted, as for the cropped deltas which causes the predicted Cp val- 
ues (zero percent leading -edge suction) to be much higher than those experimentally 
measured. The Cp data do follow Cp) q + Cl tan a (zero percent leading -e(^e suc- 
tion) where Cl is the experimental value. The Cm data are well predicted espe- 
cially at = -40°. For = 40° the extra induced effect to lift is lost, as in refer- 
ence 19 for the arrow wing, bringing the predictions of the Cl variations with a and 
Cp) into good agreement with the data. The Cm data are well predicted for lift 
coefficients up to ~0.5, after which the data are more positive. 

Flow pictures of the leading-edge shed vortex showed that breakdown occurs ahead 
of the trailing edge for the cropped -diamond, cropped -delta, and cropped-arrow wings at 
Q! ~ 25°, 24°, and 22°, respectively. The Cl variation with o; shows for the respec- 
tive models that, above these values of a, Cl increases at a slower rate than before. 

Arrow wing. - Figure 21 presents static longitudinal aerodynamic characteristics 
plotted against Cl f°i’ ^ cropped-arrow wing. In general, the data are very well pre- 
dicted up to Cl = 0.7, after which Cl falls off slightly and from the Cm data it 
appears that some or all of the contribution from the side edge is being lost. This also 
shows up in the Cl variation with a. Around Cl = 1.0 the leading -edge -vortex 
breakdown occurs ahead of the trailing edge resulting in additional loss in Cl. 

COMPARISONS WITH OTHER THEORIES 
Other Theories 

General. - Theoretical estimates for the effects of separated flow along streamwise 
side edges on the aerodynamic characteristics of swept and unswept wings have been 
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made by various authors (refs, 4 to 11). Of these the ones most commonly referenced 
are Bollay, Flax, and Gersten (refs. 4, 6, and 8, respectively). These along with the 
more recent methods of Belotserkovskii (ref. 10) and Bradley (ref. 11) will be briefly 
reviewed after which comparisons will be made with the present method and with experi- 
mental data. 

Bollay. - Bollay's flow model (ref. 4) is associated with a rectangular wing and 
assumes attached flow around the leading edge with the flow around the side edge sepa- 
rating and trailing off at some angle, nominally Of/2, to the free stream. A further 
assumption, which agrees well with experiment, is that the span loading is constant. The 
resultant formulation retains sin^ a and sin a cos a terms in the expression 

and appears adequate for the extremes of aspect ratio. However, for aspect ratios from 
1 to 6 this method overpredicts the experimental C^. Values of ACjj could be com- 
puted with Bollay’s method but, due to the employment of the mean pressure, Cm can- 
not be determined with the same accuracy. This method was not intended to be used for 
wings with sharp leading edges. 

Flax . - Flax and Lawrence (ref, 6) use a combination of attached- and cross -flow 
methods to arrive at procedures for calculating and Cm in terms of the sum of 

Q! -contributions up through for delta and rectangular wings. The coefficient of 
in the Cl equation is taken to be the two-dimensional cross -flow drag coefficient. 
Its value is 2,0 for rectangular wings with separated flow around the side edges and 
attached flow around the leading edges. The same value is used for thin delta wings. 

For rectangular wings with sharp leading and side edges the value of the two-dimensional 
cross -flow drag coefficient is replaced by an empirically determined value of 3,0. Flax 
and Lawrence state that retaining the trigonometric term in the term to make it 

sin^ a, "while more rigorously correct . , . does not always lead to results in better 
agreement with experiment," The Cm equation, which is empirically developed, has 
one linear term in a and two terms. The second term is added to contrib- 
ute a moment from a force located at the half root chord, which is found in the present 
paper to be correct only for Al = 0, Values of ACq could be computed with this 
model. 

Gersten . - Gersten (ref. 8) employs Bollay's flow model in spanwise strips which 
allows him to represent a separated flow about the leading edge for rectangular and 
delta -t 3 q)e wings. The flow trails off the surface at half the wing angle of attack result- 
ing in a full wake flow over the wing surface. Unfortunately, this does not represent well 
the physical flow model of wings with both sharp leading and side edges at least in the 
low -to -moderate angle -of -attack range where flow reattachment occurs, Gersten, like 
Flax and Lawrence, has equations for Cl and Cm which depend on a and 
(trigonometric terms neglected) where the coefficients are determined from two 


20 


contributions: (1) linear theory (attached flow) and (2) nonlinear theory (shed vorticity 
all across the wing surface). Values of ACd can be calculated easily with this method. 

Garner's method (ref. 9) was developed to be an extension of Gersten's but this 
effort was concerned primarily with establishing a similar technique having a different 
attached -flow solution. 

Belotserkovskii . - Belotserkovskii (ref. 10) solves the separated side -edge problem 
by using a discrete vortex lattice which has no trailing leg at the tip but has the bound 
filaments trail across the tip and wrap up into a shed -vortex sheet. The effect of the 
strength and location of this rolled-up tip sheet is also taken into account, by iteration, 
in obtaining a solution. The leading-edge flow is assumed to be attached. Normal -force 
and pitching -moment results are computed. 

Bradley . - Bradley's method (ref. 11) is based on computing the leading- and side- 
edge -suction forces and then employing the generalized suction analogy - similar to the 
present method^ - to predict Cl and ACd for a variety of planforms. The major 
difference between the two methods is the potential-flow solutions employed to arrive at 
the suction forces. The present method uses a modified Multhopp approach with contin- 
uous pressure loadings, whereas Bradley's method uses the discrete loadings of the 
vortex -lattice procedure. Limited comparisons made with the two methods for wings 
having separated flows indicate that they are in excellent agreement. Pitching moments 
could be determined with this method but none were given in reference 11. 

Comparisons 

Comparisons of predictions - by the previously discussed theories and the present 
method - and experimental data are to be made for representative or the complete 
series of planforms which have been analyzed and results published. 

Figure 22 presents results for = 0.2 and 1.0 rectangular wings at M « 0 from 
the present method and experiment and from the methods of Bollay and Belotserkovskii 
(jR = 1.0 only). The data for both wings are seen to be predicted better with the pres- 
ent method over the entire o? -range than with the other methods. 

Figure 23 shows static longitudinal aerodynamic characteristics plotted against a 
at M = 0 for an .R = 1.0 rectangular wing obtained from the present method and 
experiment along with those from Flax, Gersten, Garner, and Belotserkovskii. The 
Cl data are best predicted by the present method over the a -range, however, the 
Cm,le data are only best predicted by the present method up to a ~ 16°. Around this 
value of a it appears that the leading-edge vortex core has moved enough to begin to 

O 

Bradley references the original application of the present method to rectangular 

wings. 
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contribute a noticeable nose -down moment. In addition to this change it would be 
expected that the pitching moment from the side edge would be affected by the change in 
flow conditions. The present method predicted no pitching -moment contribution from 
the leading-edge vortex lift since the reference point was at the leading edge. 

Figure 24 presents the Cj^ and Cm variations with a. for two wings at M = 0 
which have A = 45°, \ = 1.0, and aspect ratios of 1 and 2. They were obtained from the 
present method and those of Gersten, Garner, and Belotserkovskii along with data from 
reference 8. At Al = 1 the present method with separated flow assumed at the leading 
and side edges predicts the Cl and Cm data more accurately than either of the other 
methods separately. For the At = 2 wing no one method works well over the entire 
0 ! -range, although up to a ~ 12° the Cl and Cm data are slightly better predicted 
by the present and Gersten's methods. A comparison of Cl p + Cl ygg with wing data 
would indicate that the leading -edge vortex lift begins to fall off noticeably for a ~ 12°. 
This could occur for several reasons, the most likely is that the shed-vortex system is 
beginning to break down near or ahead of the wing trailing edge. 

Zero Aspect-Ratio Limit 

At AR = 0.2 the present method yields a value of Ky^se =2.8 which results in 
a good estimation of Cl and Cj^ as a function of a. However, a comparison with 
Cjq- data found in reference 4 for Al < 0.2 shows that the method of Bollay (ref. 4) 
estimates more closely the experimental results and both Bollay' s estimates and the 
experimental data are smaller than the present -method predictions. Evidently, for some 
aspect ratio between 0.134 and 0.20 the flow no longer reattaches on the upper surface 
over any significant angle -of -attack range, as required by the present method, which 
results in a pure wake flow and a sharp reduction in the actual value of Ky ge- (Note 
that if Ky^se = 2.0 for a At = 0 wing a pure wake flow could be implied since this 
would be the same value as that for the cross -flow drag coefficient at this aspect ratio.) 
Gersten's results (ref. 8) indicate a nonlinear coefficient of =5.9 at Al = 0 which is 
almost twice the tt - result from the present method. 

CONCLUSIONS 

An analysis of an extension of the leading -edge -suction analogy of Polhamus (NASA 
TN D-3767) to wings with separated flow around their side edges is presented. Compar- 
isons of predicted and experimental data and other analytical methods have yielded the 
following conclusions: 

(1) Based on slender -wing theory the value of the term associated with the lift from 
the separated flow at the side edges K is tt at zero aspect ratio. 
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is zero. 


(2) The infinite -aspect -ratio value of Ky^se 

(3) For rectangular wings with aspect ratio less than 2 the value of Ky,se is 

larger than the value of the term associated with the lift at the leading edge Ky^ge- 
Above an aspect ratio of 2, Ky^ie is larger than Ky^se ior all aspect ratios below 
2,8 the totals of Ky^ge Ky ig are very near tt. 

(4) For cropped-delta wings the totals of Ky gg and Ky ig for a wide range of 

sweep angle, taper, and aspect ratios varied less than ±10 percent about n; however, both 
*^v,se ^v,le varied widely. 

(5) For cropped wings with leading- and trailing-edge sweep the total of Ky gg 
and Ky ig increases with increasing aspect ratio, leading -edge sweep, trailing-edge 
sweep (becoming more arrow like), and tip chord. 

(6) The lift variation was well predicted over an angle -of -attack range for four 
rectangular wings of aspect ratios 0.5, 1.0, 1.5, and 2.0 having round leading edges and 
vertical side edges using Ky^se as the only separated -flow lift contributor. 

(7) Experimental data of pitching moment, lift, and drag for rectangular, cropped- 
delta, cropped -diamond, and cropped-arrow wings were predicted reasonably well by the 
present method. Some drag correlations confirmed predicted drag values below those 
of 100 -percent leading -edge suction. 

(8) Without the side -edge contribution to separated -flow lift and moment, the value 
of predicted lift would, in general, be too low and the moments could be more positive. 

(9) Sweptback wings which have large amounts of area behind the point of maximum 
span develop lift values in excess of those predicted due to additional induced effects 
associated with the actual shed -vortex system. 

(10) It has been determined that the present method predicts wing static longitudinal 
aerodynamic data as well as, or better than, many other methods, 

(11) The separated -flow model assumed by the present method is more nearly con- 
sistent with the real flow than that of some other methods. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 23, 1974. 
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TABLE I. - MODEL DIMENSIONS 


Rectangular wings 

JR 

b 

Cp 3.nd, ^x*6f 

Sref 

xref 

cm 

in. 

cm 

in. 

m'^ 

ft^ 

cm 

in. 

0.20 

.30 

.40 

1.00 

3.00 

25.40 

30.48 

30.48 

50.80 

76.20 

10 

12 

12 

20 

30 

127.00 

101.60 

76.20 

50.80 

25.40 

50 

40 

30 

20 

10 

0.322 

.309 

.232 

.258 

.193 

3.47 

3.33 

2.50 

2.78 

2.08 

31.75 

25.40 

19.05 

12.70 

6.35 

12.5 

10.0 

7.5 
5.0 

2.5 

Notched rectangular wing 


b 


°t 

Cref 

bn 

Sref 

Xref 


cm 

in. 

cm 

in. 

cm 

in. 

cm 

in. 

B 

in. 

m2 

r2 


m 

1.00 

50.80 

20 

38.94 

15.33 


m 

51.99 

20.47 

15.24 

6 

0.258 

m 

13.00 

m 

Cropped-delta wings with Cr = 25.4 cm (10 in.) and A = 63° (tested at M = 0.6) 

(see ref. 18) 

JR 


b 

Cref 

Sref 

Xref 

cm 



■a 



m 

in. 

1.64 

1.33 

1.08 

.86 

1 

22.89 

20.32 

17.80 

15.27 

9.01 
8.00 

7.01 

6.01 


6.73 

6.89 

7.13 

7.43 



12.57 

12.27 

11.81 

11.25 

4.95 

4.83 

4.65 

4.43 

Cropped wings with ct = 33.22 cm (13.08 in.) and A = 63° (tested at M = 0.2) 

JR 

deg 

b 

Cp 

Cref 

Sref 

Xref 

cm 

in. 

cm 

in. 

cm 

in. 

m2 

ft2 

cm 

in. 

0.87 

.74 

1.07 


■ 

I 

83.08 

104.39 

62.31 

32.71 

41.10 

24.53 

61.72 

74.96 

48.92 

24.30 

29.51 

19.26 

0.295 

.349 

.242 

3.18 

3.76 

2.60 

36.80 

39.37 

34.67 

14.49 

15.50 
13.65 


Ct = 11.81 cm (4.65 in.) and A = 63° (tested at M = 0.2) 

JR 

D, 

deg 

b 


Cref 

Sref 

Xref 

■i 


cm 

in. 

cm 

in. 


ft2 


n 

1.92 

38.86 

50.80 

20 

41.20 

16.22 

29.21 

11.50 

0.135 

1.45 

27.61 

10.87 
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Figure 2.- Effect of chord and span loadings on K of rectangular wings 
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Figure 4.- Effect of aspect ratio on tip -suction distribution of rectangular wings at 





Figure 5.- Variation of jSK , 


/3K^ le’ tot aspect ratio and Mach number 

for rectangular wings. 
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Figure 7.- Effect of leading -edge sweepback and Mach number on Ky^le and Ky ge for cropped-delta wings. 






A = 63® 



Figure 9.- Effect of leading -edge sweepback and tip chord on Ky^tot selected 

families of cropped wings at M = 0. 


35 










Photographs of the sharp -edged cropped wing models tested. 
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Figure 11.- Continued. 
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Figure 11.- Continued, 



(d) Cropped arrow (2), 
Figure 11.- Concluded. 
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(a) Three-quarter rear view. 

Figure 12.- Typical wind-tunnel setup for a sharp-edged unswept wing. 



Figure 12.- Concluded 
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(a) Three-quarter rear view. 

Figure 13.- Typical wind-tunnel setup for a sharp-edged cropped wing. 
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Figure 15.- Effect of aspect ratio on Cj^ for rectangular wings with round leading 
edges and vertical sides at M ~ 0 (Gottingen 409 airfoil; thickness ratio = 0.126) 
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Figure 16.- Drag variation of rectangular wing associated with different edge condition assumptions at 



Flow 







51 


Figure 16.- Continued. 
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Figure 16.- Concluded. 
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(a) M = 0.2. 

Figure 17.- Effect of aspect ratio on some static longitudinal aerodynamic 
characteristics for sharp -edged rectangular wings at M = 0.2. 
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(e) m = 3.0. 
Figure 17.- Concluded. 
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Figure 18.- Some static longitudinal aerodynamic characteristics of a sharp-edged 
notched rectangular wing. JR = 1.0; M = 0.2; X = 1.435. 



Cl 

(a) X = 0.10. 

Figure 19,- Effect of taper ratio on some static longitudinal aerodynamic characteristics 
for cropped delta wings with A = 63° and NACA 63A002 airfoils at M = 0.6. Experi- 
ment obtained from reference 18. 
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Cl 

(b) \ = 0.20. 
Figure 19.- Continued. 
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(d) \ = 0.40. 
Figure 19.- Concluded. 




(a) n = -40°. 

Figure 20.- Effect of trailing -edge sweep on some static longitudinal aerodynamic 
characteristics for sharp -edged wings with constant tip chord \o-o = 0.4, and 
A = 630 at m = 0.2. 
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(b) a = 0°. 


Figure 20,- Continued 
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Figure 20,- Concluded. 
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Figure 23.- Theoretical and experimental static longitudinal aerodynamic char 
acteristics for a sharp edged rectangular wing of ^ = 1.0 at M ~ 0. 
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Figure 24.- Effect of aspect ratio on the theoretical and experimental static longitudinal 
aerod5mamic characteristics for two wings. A = 45°; A. = 1.0; M = 0. 
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